
R : Copyright 2003, The R Foundation for Statistical Computing
Version 1.8.1  (2003-11-21), ISBN 3-900051-00-3

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for a HTML browser interface to help.
Type 'q()' to quit R.

> invisible(options(echo = TRUE))
> #
> # This estimates the Florida model for five categories
> 
> options(width=120);
> set.seed(19342)
> 
> if (!exists("multinomRob"))
+ 	library(multinomRob)
Loading required package: rgenoud 
Loading required package: MASS 
Loading required package: mvtnorm 
> 
> ## read data
> dta <- read.table("REPLICATION_DATA.asc", header=TRUE)
> 
> dta$pdemnader96 <- dta$pdem96 + dta$pnader96 ;
> dta$other00 <- dta$other.candidates00 ;
> 
> ## See names
> names(dta)
 [1] "name"               "TotalVotes00"       "TotalVotes96"       "pDemRV00"           "pRepRV00"          
 [6] "pBlackDemRV00"      "pNonBlackDemRV00"   "pDemRV96"           "pRepRV96"           "pChangeDemRV00"    
[11] "pChangeRepRV00"     "pdem00"             "prep00"             "buchanan00"         "nader00"           
[16] "gore00"             "bush00"             "other.candidates00" "pperot96"           "pnader96"          
[21] "pdem96"             "prep96"             "pGreenRV00"         "pChangeDemRV00.1"   "pChangeRepRV00.1"  
[26] "pcuban00"           "pr1...1."           "pr1...2."           "pr1...3."           "pr1...4."          
[31] "pr2...1."           "pr2...2."           "pr2...3."           "pr2...4."           "pr3...1."          
[36] "pr3...2."           "pr3...3."           "pr3...4."           "inc"                "cen00.native"      
[41] "cen00.asian"        "college"            "cen00.black"        "cen00.hispanic"     "cen00.white"       
[46] "density"            "cen00.population"   "pdemnader96"        "other00"           
> ## And define a model using these names:
> model <- list(buchanan00 ~ pperot96 + pChangeRepRV00 + pcuban00 + pr1...1.,
+               nader00    ~ pdemnader96 +  pGreenRV00 + pcuban00 +   pr1...2.,
+               gore00     ~ pdem96 + pChangeDemRV00 + pcuban00 + pr1...3.,
+               bush00     ~ prep96 + pChangeRepRV00 + pcuban00 + pr1...4.,
+               other00    ~ 0
+               )
> 
> ## Set parameters for Genoud
> zz.genoud.parms <- list(wait.generations      = 30,
+                         pop.size              = 2000,
+                         P5                    = 500,
+                         scale.domains         = 4
+                         )
> 
> load("FLfull1.RData");
> starts <- muFLfull1$genoud$par;
> 
> muFLfull1 <- multinomRob(model,
+                  data=dta,
+                  starting.values = starts,
+                  genoud.parms   = zz.genoud.parms,
+                  print.level    = 3 ) 

Your Model (xvec):
                                                           buchanan00 nader00 gore00 bush00 other00
(Intercept)/(Intercept)/(Intercept)/(Intercept)/NA                  1       1      1      1       0
pperot96/pdemnader96/pdem96/prep96/NA                               1       1      1      1       0
pChangeRepRV00/pGreenRV00/pChangeDemRV00/pChangeRepRV00/NA          1       1      1      1       0
pcuban00/pcuban00/pcuban00/pcuban00/NA                              1       1      1      1       0
pr1...1./pr1...2./pr1...3./pr1...4./NA                              1       1      1      1       0
multinomRob(): Starting Values 
                                                            buchanan00      nader00     gore00      bush00 other00
(Intercept)/(Intercept)/(Intercept)/(Intercept)/NA         -0.60538069   0.76560521 3.33351379  4.09429714       0
pperot96/pdemnader96/pdem96/prep96/NA                       4.64215952   1.14520548 3.07427094  1.70258490       0
pChangeRepRV00/pGreenRV00/pChangeDemRV00/pChangeRepRV00/NA 13.06463688 660.40257319 1.01749984  2.64429532       0
pcuban00/pcuban00/pcuban00/pcuban00/NA                     -0.88909869  -1.05530975 1.46343604  2.76939817       0
pr1...1./pr1...2./pr1...3./pr1...4./NA                     -0.04636576   0.04190072 0.01974885 -0.01148855       0
multinomRob(): starting fit = 5.046338 


Fri Mar  5 00:29:17 2004
Domains:
 -3.026903e+00   <=  X1   <=    1.816142e+00 
 -1.392648e+01   <=  X2   <=    2.321080e+01 
 -3.919391e+01   <=  X3   <=    6.532318e+01 
 -4.445493e+00   <=  X4   <=    2.667296e+00 
 -2.318288e-01   <=  X5   <=    1.390973e-01 
 -2.296816e+00   <=  X6   <=    3.828026e+00 
 -3.435616e+00   <=  X7   <=    5.726027e+00 
 -1.981208e+03   <=  X8   <=    3.302013e+03 
 -5.276549e+00   <=  X9   <=    3.165929e+00 
 -1.257022e-01   <=  X10  <=    2.095036e-01 
 -1.000054e+01   <=  X11  <=    1.666757e+01 
 -9.222813e+00   <=  X12  <=    1.537135e+01 
 -3.052500e+00   <=  X13  <=    5.087499e+00 
 -4.390308e+00   <=  X14  <=    7.317180e+00 
 -5.924656e-02   <=  X15  <=    9.874427e-02 
 -1.228289e+01   <=  X16  <=    2.047149e+01 
 -5.107755e+00   <=  X17  <=    8.512925e+00 
 -7.932886e+00   <=  X18  <=    1.322148e+01 
 -8.308195e+00   <=  X19  <=    1.384699e+01 
 -5.744275e-02   <=  X20  <=    3.446565e-02 

Data Type: Floating Point
Operators (code number, name, population) 
	(1) Cloning........................... 	116
	(2) Uniform Mutation.................. 	118
	(3) Boundary Mutation................. 	118
	(4) Non-Uniform Mutation.............. 	118
	(5) Polytope Crossover................ 	1175
	(6) Multiple Point Simple Crossover... 	118
	(7) Whole Non-Uniform Mutation........ 	118
	(8) Heuristic Crossover............... 	118
	(9) Local-Minimum Crossover........... 	0

SOFT Maximum Number of Generations: 100
Maximum Nonchanging Generations: 30
Population size       : 2000
Convergence Tolerance: 1.000000e-07

Using the BFGS Derivative Based Optimizer on the Best Individual Each Generation.
Not Checking Gradients before Stoping
Using Out of Bounds Individuals.



Minimization Problem.



STARTING VALUES


The 2 best initial individuals are
-0.61	 4.64	13.06	-0.89	-0.05	 0.77	 1.15	660.40	-1.06	 0.04	 3.33	 3.07	 1.02	 1.46	 0.02	 4.09	 1.70	 2.64	 2.77	-0.01	
fitness = 5.046338e+00

-0.38	-8.13	 7.43	-1.73	-0.13	 1.16	-1.61	-606.95	-2.05	-0.09	-8.66	 2.79	 0.41	-3.16	 0.07	 3.80	 2.03	 7.51	-5.30	-0.01	
fitness = 1.868610e+01


The worst fit of the population is: 6.593387e+06




Generation#	    Solution Value

      0 	5.046338e+00
      1 	5.046338e+00
      2 	5.046338e+00
      3 	5.046338e+00
      4 	5.046338e+00
      5 	5.046338e+00
      6 	5.046338e+00
      7 	5.046338e+00
      8 	5.046338e+00
      9 	5.046338e+00
     10 	5.046338e+00
     11 	5.046338e+00
     12 	5.046338e+00
     13 	5.046338e+00
     14 	5.046338e+00
     15 	5.046338e+00
     16 	5.046338e+00
     17 	5.046338e+00
     18 	5.046338e+00
     19 	5.046338e+00
     20 	5.046338e+00
     21 	5.046338e+00
     22 	5.046338e+00
     23 	5.046338e+00
     24 	5.046338e+00
     25 	5.046338e+00
     26 	5.046338e+00
     27 	5.046338e+00
     28 	5.046338e+00
     29 	5.046337e+00
     30 	5.046336e+00
     31 	5.046335e+00
     32 	5.046335e+00
     33 	5.046335e+00
     34 	5.046335e+00
     35 	5.046335e+00
     36 	5.046335e+00
     37 	5.046335e+00
     38 	5.046335e+00
     39 	5.046335e+00
     40 	5.046335e+00
     41 	5.046335e+00
     42 	5.046335e+00
     43 	5.046335e+00
     44 	5.046335e+00
     45 	5.046335e+00
     46 	5.046335e+00
     47 	5.046335e+00
     48 	5.046335e+00
     49 	5.046335e+00
     50 	5.046335e+00
     51 	5.046335e+00
     52 	5.046335e+00
     53 	5.046335e+00
     54 	5.046335e+00
     55 	5.046335e+00
     56 	5.046333e+00
     57 	5.046324e+00
     58 	5.046323e+00
     59 	5.046318e+00
     60 	5.046318e+00
     61 	5.046317e+00
     62 	5.046315e+00
     63 	5.046315e+00
     64 	5.046315e+00
     65 	5.046315e+00
     66 	5.046315e+00
     67 	5.046315e+00
     68 	5.046315e+00
     69 	5.046313e+00
     70 	5.046310e+00
     71 	5.046309e+00
     72 	5.046308e+00
     73 	5.046306e+00
     74 	5.046305e+00
     75 	5.046304e+00
     76 	5.046304e+00
     77 	5.046304e+00
     78 	5.046304e+00
     79 	5.046304e+00
XXX
XXX
XXX
XXX
XXX
XXX
     80 	5.046301e+00
     81 	5.046301e+00
     82 	5.046300e+00
     83 	5.046300e+00
     84 	5.046300e+00
     85 	5.046300e+00
     86 	5.046300e+00
     87 	5.046300e+00
     88 	5.046300e+00
     89 	5.046300e+00
     90 	5.046300e+00
     91 	5.046300e+00
     92 	5.046300e+00
     93 	5.046300e+00
     94 	5.046300e+00
     95 	5.046300e+00
     96 	5.046300e+00
     97 	5.046300e+00
     98 	5.046300e+00
     99 	5.046300e+00
    100 	5.046300e+00

Increasing Soft Maximum Generation Limit by 100 (MaxGenerations) to 200.
I'm doing this because the fitness is still improving.

    101 	5.046300e+00
    102 	5.046300e+00
    103 	5.046300e+00
    104 	5.046300e+00
    105 	5.046300e+00
    106 	5.046300e+00
    107 	5.046300e+00
    108 	5.046300e+00
    109 	5.046300e+00
    110 	5.046300e+00
    111 	5.046300e+00
    112 	5.046300e+00
    113 	5.046300e+00
    114 	5.046300e+00
    115 	5.046300e+00
    116 	5.046300e+00
    117 	5.046300e+00
    118 	5.046300e+00
    119 	5.046300e+00
    120 	5.046300e+00
    121 	5.046300e+00
    122 	5.046300e+00
    123 	5.046300e+00
    124 	5.046300e+00
    125 	5.046300e+00
    128 	5.046300e+00
    134 	5.046300e+00
    135 	5.046300e+00
    136 	5.046300e+00
    137 	5.046300e+00
    138 	5.046300e+00

Soft Generation Wait Limit Hit.
No Improvement in 30 Generations

Best Fit Found at Generation 138
Fit Value = 5.046300e+00


Parameters at the Solution (value, gradient):

 X[ 1] :	-6.053806e-01	G[ 1] :	3.457535e+00
 X[ 2] :	4.642160e+00	G[ 2] :	-4.471773e-01
 X[ 3] :	1.306465e+01	G[ 3] :	-1.888999e+00
 X[ 4] :	-8.891408e-01	G[ 4] :	2.722228e-01
 X[ 5] :	-4.636573e-02	G[ 5] :	-6.751346e+01
 X[ 6] :	7.656044e-01	G[ 6] :	-2.444179e-01
 X[ 7] :	1.145206e+00	G[ 7] :	-1.065650e+01
 X[ 8] :	6.604027e+02	G[ 8] :	-1.461769e-04
 X[ 9] :	-1.055294e+00	G[ 9] :	-6.399390e-03
 X[10] :	4.190320e-02	G[10] :	-1.640088e+00
 X[11] :	3.333514e+00	G[11] :	-1.823503e+02
 X[12] :	3.074271e+00	G[12] :	-7.510604e+01
 X[13] :	1.017515e+00	G[13] :	3.377232e+00
 X[14] :	1.463459e+00	G[14] :	-1.454379e-01
 X[15] :	1.974802e-02	G[15] :	-1.254962e+01
 X[16] :	4.094297e+00	G[16] :	1.237627e+02
 X[17] :	1.702585e+00	G[17] :	6.346459e+01
 X[18] :	2.644293e+00	G[18] :	-2.723467e-01
 X[19] :	2.769398e+00	G[19] :	-8.255148e-02
 X[20] :	-1.148819e-02	G[20] :	-1.413830e+01

Solution Found Generation 138
Number of Generations Run 150


Fri Mar  5 00:36:40 2004
Total run time : 0 hours 7 minutes and 23 seconds


LQD Results:
                                                            buchanan00     nader00     gore00      bush00 other00
(Intercept)/(Intercept)/(Intercept)/(Intercept)/NA         -0.60538064   0.7656044 3.33351369  4.09429719       0
pperot96/pdemnader96/pdem96/prep96/NA                       4.64215952   1.1452056 3.07427089  1.70258489       0
pChangeRepRV00/pGreenRV00/pChangeDemRV00/pChangeRepRV00/NA 13.06464774 660.4027146 1.01751505  2.64429298       0
pcuban00/pcuban00/pcuban00/pcuban00/NA                     -0.88914078  -1.0552943 1.46345868  2.76939809       0
pr1...1./pr1...2./pr1...3./pr1...4./NA                     -0.04636573   0.0419032 0.01974802 -0.01148819       0

LQD sigma: 5.0463 


(multinomTanh):
Tanh Estimates
                                                            buchanan00      nader00     gore00     bush00 other00
(Intercept)/(Intercept)/(Intercept)/(Intercept)/NA         -0.60456197   0.73715558 3.34631760  4.0883995       0
pperot96/pdemnader96/pdem96/prep96/NA                       4.62711626   1.23761373 3.05828092  1.7197376       0
pChangeRepRV00/pGreenRV00/pChangeDemRV00/pChangeRepRV00/NA 13.06013844 644.39215934 1.05182337  2.6300305       0
pcuban00/pcuban00/pcuban00/pcuban00/NA                     -1.06089338  -0.97996697 1.47728543  2.7618678       0
pr1...1./pr1...2./pr1...3./pr1...4./NA                     -0.05097921   0.03514663 0.02062767 -0.0086593       0

Tanh Sandwich SEs
                                                           buchanan00     nader00     gore00     bush00 other00
(Intercept)/(Intercept)/(Intercept)/(Intercept)/NA         0.14605079  0.18273414 0.15279549 0.14074436       0
pperot96/pdemnader96/pdem96/prep96/NA                      1.31181220  0.45275434 0.31574520 0.29157093       0
pChangeRepRV00/pGreenRV00/pChangeDemRV00/pChangeRepRV00/NA 1.23328550 77.47034914 0.92423871 0.71561997       0
pcuban00/pcuban00/pcuban00/pcuban00/NA                     0.45294960  0.27248483 0.24199806 0.22598938       0
pr1...1./pr1...2./pr1...3./pr1...4./NA                     0.02428515  0.01437326 0.01043674 0.00747727       0

TANH sigma: 4.349402 


**********************************

Iteration: 1 
using best non-genoud LQD sigma. Provided by Tanh.


Fri Mar  5 00:36:41 2004
Domains:
 -1.813686e+00   <=  X1   <=    6.045620e-01 
 -4.627116e+00   <=  X2   <=    1.388135e+01 
 -1.306014e+01   <=  X3   <=    3.918042e+01 
 -3.182680e+00   <=  X4   <=    1.060893e+00 
 -1.529376e-01   <=  X5   <=    5.097921e-02 
 -7.371556e-01   <=  X6   <=    2.211467e+00 
 -1.237614e+00   <=  X7   <=    3.712841e+00 
 -6.443922e+02   <=  X8   <=    1.933176e+03 
 -2.939901e+00   <=  X9   <=    9.799670e-01 
 -3.514663e-02   <=  X10  <=    1.054399e-01 
 -3.346318e+00   <=  X11  <=    1.003895e+01 
 -3.058281e+00   <=  X12  <=    9.174843e+00 
 -1.051823e+00   <=  X13  <=    3.155470e+00 
 -1.477285e+00   <=  X14  <=    4.431856e+00 
 -2.062767e-02   <=  X15  <=    6.188302e-02 
 -4.088400e+00   <=  X16  <=    1.226520e+01 
 -1.719738e+00   <=  X17  <=    5.159213e+00 
 -2.630030e+00   <=  X18  <=    7.890091e+00 
 -2.761868e+00   <=  X19  <=    8.285603e+00 
 -2.597790e-02   <=  X20  <=    8.659300e-03 

Data Type: Floating Point
Operators (code number, name, population) 
	(1) Cloning........................... 	116
	(2) Uniform Mutation.................. 	118
	(3) Boundary Mutation................. 	118
	(4) Non-Uniform Mutation.............. 	118
	(5) Polytope Crossover................ 	1175
	(6) Multiple Point Simple Crossover... 	118
	(7) Whole Non-Uniform Mutation........ 	118
	(8) Heuristic Crossover............... 	118
	(9) Local-Minimum Crossover........... 	0

SOFT Maximum Number of Generations: 100
Maximum Nonchanging Generations: 30
Population size       : 2000
Convergence Tolerance: 1.000000e-07

Using the BFGS Derivative Based Optimizer on the Best Individual Each Generation.
Not Checking Gradients before Stoping
Using Out of Bounds Individuals.



Minimization Problem.



STARTING VALUES


The 2 best initial individuals are
-0.60	 4.63	13.06	-1.06	-0.05	 0.74	 1.24	644.39	-0.98	 0.04	 3.35	 3.06	 1.05	 1.48	 0.02	 4.09	 1.72	 2.63	 2.76	-0.01	
fitness = 5.216629e+00

-0.76	-0.07	10.35	-3.05	-0.10	 0.02	 2.77	-4.76	-2.76	-0.02	 0.45	 8.42	 0.25	 1.77	 0.03	 4.37	 1.02	 2.92	-2.29	-0.01	
fitness = 1.225279e+01


The worst fit of the population is: 7.725033e+04




Generation#	    Solution Value

      0 	5.216629e+00
      1 	5.109284e+00
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
      2 	5.093736e+00
      3 	5.082081e+00
XXX
      4 	5.075538e+00
XXX
      5 	5.063506e+00
      6 	5.061159e+00
      7 	5.059186e+00
      8 	5.058716e+00
      9 	5.057262e+00
     10 	5.054946e+00
     11 	5.050150e+00
XXX
XXX
     12 	5.041680e+00
     13 	5.041420e+00
XXX
XXX
     14 	5.041106e+00
     15 	5.021111e+00
     16 	5.021092e+00
     17 	5.020992e+00
     18 	5.019785e+00
     19 	5.019758e+00
     20 	5.019148e+00
     21 	5.011730e+00
     22 	5.003014e+00
     23 	4.988747e+00
     24 	4.986649e+00
     25 	4.983291e+00
     26 	4.978127e+00
     27 	4.972030e+00
     28 	4.971777e+00
     29 	4.969598e+00
     30 	4.969274e+00
     31 	4.968876e+00
     33 	4.968696e+00
     34 	4.968492e+00
     35 	4.968201e+00
     36 	4.968070e+00
     37 	4.967830e+00
     38 	4.967810e+00
     39 	4.967801e+00
XXX
XXX
XXX
     40 	4.967781e+00
XXX
XXX
     41 	4.967455e+00
XXX
XXX
XXX
XXX
     43 	4.967414e+00
     44 	4.967410e+00
     45 	4.967409e+00
     46 	4.967409e+00
     47 	4.967409e+00
     48 	4.967408e+00
     49 	4.967407e+00
     50 	4.967405e+00
     51 	4.967399e+00
     52 	4.967398e+00
     53 	4.967395e+00
     54 	4.967394e+00
     55 	4.967392e+00
     56 	4.967391e+00
     57 	4.967391e+00
     58 	4.967390e+00
     59 	4.967390e+00
     60 	4.967390e+00
     61 	4.967390e+00
     62 	4.967390e+00
     63 	4.967389e+00
     64 	4.967388e+00
     65 	4.967387e+00
     66 	4.967386e+00
     67 	4.967385e+00
     68 	4.967377e+00
     69 	4.967339e+00
     70 	4.967221e+00
     71 	4.966903e+00
     72 	4.966857e+00
XXX
XXX
XXX
     73 	4.966711e+00
     74 	4.966676e+00
     75 	4.966651e+00
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     76 	4.966529e+00
     77 	4.966465e+00
     78 	4.966432e+00
     79 	4.966425e+00
     80 	4.966414e+00
     81 	4.966410e+00
     82 	4.966396e+00
     83 	4.966392e+00
     84 	4.966390e+00
XXX
XXX
     85 	4.966377e+00
     86 	4.966376e+00
     87 	4.966376e+00
     88 	4.966374e+00
     89 	4.966374e+00
     90 	4.966374e+00
XXX
XXX
     91 	4.966373e+00
     92 	4.966373e+00
     93 	4.966373e+00
     94 	4.966373e+00
     95 	4.966373e+00
     96 	4.966373e+00
     97 	4.966372e+00
     98 	4.966372e+00
     99 	4.966372e+00
    100 	4.966372e+00

Increasing Soft Maximum Generation Limit by 100 (MaxGenerations) to 200.
I'm doing this because the fitness is still improving.

    101 	4.966372e+00
    102 	4.966372e+00
    103 	4.966371e+00
    104 	4.966371e+00
    105 	4.966371e+00
    106 	4.966371e+00
    107 	4.966371e+00
    108 	4.966370e+00
    109 	4.966370e+00
    110 	4.966370e+00
    111 	4.966370e+00
    112 	4.966369e+00
    113 	4.966369e+00
    114 	4.966368e+00
    115 	4.966368e+00
    116 	4.966368e+00
    117 	4.966368e+00
    118 	4.966368e+00
    119 	4.966368e+00
    120 	4.966368e+00
    121 	4.966368e+00
    122 	4.966368e+00
    123 	4.966367e+00
    124 	4.966367e+00
    125 	4.966367e+00
    126 	4.966367e+00
    127 	4.966367e+00
    128 	4.966367e+00
    129 	4.966367e+00
    130 	4.966367e+00
    131 	4.966367e+00
    132 	4.966367e+00
    133 	4.966367e+00
    134 	4.966367e+00
    135 	4.966367e+00
    136 	4.966367e+00
    137 	4.966367e+00
    138 	4.966367e+00
    139 	4.966367e+00
    140 	4.966367e+00
    141 	4.966367e+00
    142 	4.966367e+00
    143 	4.966367e+00
    144 	4.966367e+00
    145 	4.966367e+00
    146 	4.966367e+00
    147 	4.966367e+00
    148 	4.966367e+00
    149 	4.966367e+00
    150 	4.966367e+00
    151 	4.966367e+00
    152 	4.966367e+00
    153 	4.966367e+00
    154 	4.966367e+00
    155 	4.966367e+00
    156 	4.966367e+00
    157 	4.966367e+00
    158 	4.966367e+00
    159 	4.966367e+00
    160 	4.966367e+00

Soft Generation Wait Limit Hit.
No Improvement in 30 Generations

Best Fit Found at Generation 160
Fit Value = 4.966367e+00


Parameters at the Solution (value, gradient):

 X[ 1] :	-6.075802e-01	G[ 1] :	1.409863e+00
 X[ 2] :	4.946807e+00	G[ 2] :	-1.347928e+00
 X[ 3] :	1.377232e+01	G[ 3] :	-2.007702e-01
 X[ 4] :	-8.358977e-01	G[ 4] :	6.041400e+00
 X[ 5] :	-4.735228e-02	G[ 5] :	-7.279250e+00
 X[ 6] :	7.454592e-01	G[ 6] :	-1.647200e+00
 X[ 7] :	1.238929e+00	G[ 7] :	-2.680514e+01
 X[ 8] :	6.443892e+02	G[ 8] :	2.260400e-05
 X[ 9] :	-1.117820e+00	G[ 9] :	-1.198830e+00
 X[10] :	4.413117e-02	G[10] :	-1.176806e+02
 X[11] :	3.344684e+00	G[11] :	2.650253e+00
 X[12] :	3.057769e+00	G[12] :	-2.595503e+01
 X[13] :	1.078261e+00	G[13] :	3.987473e+00
 X[14] :	1.473881e+00	G[14] :	-7.750409e+01
 X[15] :	2.304997e-02	G[15] :	-1.387442e+03
 X[16] :	4.088322e+00	G[16] :	2.904776e+00
 X[17] :	1.721391e+00	G[17] :	4.775448e+00
 X[18] :	2.629912e+00	G[18] :	5.904195e-01
 X[19] :	2.795618e+00	G[19] :	3.180017e+01
 X[20] :	-9.560690e-03	G[20] :	5.266799e+01

Solution Found Generation 160
Number of Generations Run 160


Fri Mar  5 00:45:46 2004
Total run time : 0 hours 9 minutes and 5 seconds

(iter 1 ): LQD Results:
                                                            buchanan00      nader00     gore00      bush00 other00
(Intercept)/(Intercept)/(Intercept)/(Intercept)/NA         -0.60758019   0.74545919 3.34468426  4.08832219       0
pperot96/pdemnader96/pdem96/prep96/NA                       4.94680671   1.23892882 3.05776912  1.72139117       0
pChangeRepRV00/pGreenRV00/pChangeDemRV00/pChangeRepRV00/NA 13.77231667 644.38920836 1.07826062  2.62991218       0
pcuban00/pcuban00/pcuban00/pcuban00/NA                     -0.83589773  -1.11782027 1.47388120  2.79561780       0
pr1...1./pr1...2./pr1...3./pr1...4./NA                     -0.04735228   0.04413117 0.02304997 -0.00956069       0

(iter 1 ): LQD sigma: 4.966367 

(iter 1 ): Tanh Estimates
                                                            buchanan00      nader00     gore00       bush00 other00
(Intercept)/(Intercept)/(Intercept)/(Intercept)/NA         -0.60289980   0.72474965 3.35097285  4.082978517       0
pperot96/pdemnader96/pdem96/prep96/NA                       4.67113304   1.28489046 3.04973640  1.732836788       0
pChangeRepRV00/pGreenRV00/pChangeDemRV00/pChangeRepRV00/NA 13.20404872 638.40858219 1.08822419  2.613743540       0
pcuban00/pcuban00/pcuban00/pcuban00/NA                     -1.05402855  -1.02616980 1.47292387  2.763455540       0
pr1...1./pr1...2./pr1...3./pr1...4./NA                     -0.05072537   0.03589573 0.02198428 -0.007951595       0

(iter 1 ): Tanh Sandwich SEs
                                                           buchanan00     nader00     gore00      bush00 other00
(Intercept)/(Intercept)/(Intercept)/(Intercept)/NA         0.14672920  0.18112582 0.15047397 0.138331574       0
pperot96/pdemnader96/pdem96/prep96/NA                      1.31495608  0.44984818 0.31071514 0.286459543       0
pChangeRepRV00/pGreenRV00/pChangeDemRV00/pChangeRepRV00/NA 1.23568813 76.62043174 0.91382444 0.709923940       0
pcuban00/pcuban00/pcuban00/pcuban00/NA                     0.44845600  0.27530268 0.23635670 0.220814772       0
pr1...1./pr1...2./pr1...3./pr1...4./NA                     0.02403404  0.01410324 0.01033199 0.007427573       0

(iter 1 ): TANH sigma: 4.318946 


**********************************

Iteration: 2 
using LQD sigma provided by Tanh.


Fri Mar  5 00:45:47 2004
Domains:
 -1.808699e+00   <=  X1   <=    6.028998e-01 
 -4.671133e+00   <=  X2   <=    1.401340e+01 
 -1.320405e+01   <=  X3   <=    3.961215e+01 
 -3.162086e+00   <=  X4   <=    1.054029e+00 
 -1.521761e-01   <=  X5   <=    5.072537e-02 
 -7.247497e-01   <=  X6   <=    2.174249e+00 
 -1.284890e+00   <=  X7   <=    3.854671e+00 
 -6.384086e+02   <=  X8   <=    1.915226e+03 
 -3.078509e+00   <=  X9   <=    1.026170e+00 
 -3.589573e-02   <=  X10  <=    1.076872e-01 
 -3.350973e+00   <=  X11  <=    1.005292e+01 
 -3.049736e+00   <=  X12  <=    9.149209e+00 
 -1.088224e+00   <=  X13  <=    3.264673e+00 
 -1.472924e+00   <=  X14  <=    4.418772e+00 
 -2.198428e-02   <=  X15  <=    6.595284e-02 
 -4.082979e+00   <=  X16  <=    1.224894e+01 
 -1.732837e+00   <=  X17  <=    5.198510e+00 
 -2.613744e+00   <=  X18  <=    7.841231e+00 
 -2.763456e+00   <=  X19  <=    8.290367e+00 
 -2.385479e-02   <=  X20  <=    7.951595e-03 

Data Type: Floating Point
Operators (code number, name, population) 
	(1) Cloning........................... 	116
	(2) Uniform Mutation.................. 	118
	(3) Boundary Mutation................. 	118
	(4) Non-Uniform Mutation.............. 	118
	(5) Polytope Crossover................ 	1175
	(6) Multiple Point Simple Crossover... 	118
	(7) Whole Non-Uniform Mutation........ 	118
	(8) Heuristic Crossover............... 	118
	(9) Local-Minimum Crossover........... 	0

SOFT Maximum Number of Generations: 100
Maximum Nonchanging Generations: 30
Population size       : 2000
Convergence Tolerance: 1.000000e-07

Using the BFGS Derivative Based Optimizer on the Best Individual Each Generation.
Not Checking Gradients before Stoping
Using Out of Bounds Individuals.



Minimization Problem.



STARTING VALUES


The 2 best initial individuals are
-0.60	 4.67	13.20	-1.05	-0.05	 0.72	 1.28	638.41	-1.03	 0.04	 3.35	 3.05	 1.09	 1.47	 0.02	 4.08	 1.73	 2.61	 2.76	-0.01	
fitness = 5.184236e+00

-0.70	 9.70	 5.07	-1.93	-0.08	 0.82	 1.28	251.30	-1.17	-0.03	 3.64	 2.48	-0.41	 3.94	-0.02	 5.02	 0.46	-0.93	-0.24	-0.01	
fitness = 9.669563e+00


The worst fit of the population is: 6.549789e+04




Generation#	    Solution Value

      0 	5.184236e+00
      1 	5.088288e+00
XXX
XXX
      2 	5.082660e+00
      3 	5.077875e+00
      4 	5.063585e+00
XXX
      5 	5.059428e+00
      6 	5.056512e+00
XXX
XXX
XXX
XXX
      7 	5.052272e+00
XXX
XXX
XXX
XXX
XXX
XXX
      8 	5.049978e+00
      9 	5.049550e+00
     10 	5.048910e+00
     11 	5.048301e+00
     12 	5.047757e+00
     13 	5.047373e+00
     14 	5.031654e+00
     15 	5.028126e+00
     16 	5.027927e+00
     17 	5.027440e+00
     18 	5.027205e+00
     19 	5.026988e+00
     20 	5.026789e+00
     21 	5.026752e+00
     22 	5.026749e+00
     23 	5.026504e+00
     24 	5.026230e+00
     25 	5.025918e+00
     27 	5.025839e+00
     28 	5.024988e+00
     29 	5.024925e+00
     30 	5.024925e+00
     31 	5.024920e+00
     32 	5.024722e+00
     33 	5.024650e+00
     34 	5.024595e+00
     35 	5.024453e+00
     36 	5.024399e+00
     37 	5.024290e+00
     38 	5.023401e+00
XXX
XXX
XXX
XXX
     39 	5.022886e+00
     40 	5.021700e+00
XXX
XXX
XXX
XXX
XXX
XXX
     41 	5.021495e+00
XXX
XXX
     42 	5.021473e+00
     43 	5.021260e+00
     44 	5.021211e+00
     45 	5.020916e+00
     46 	5.020229e+00
     47 	5.020220e+00
XXX
XXX
XXX
XXX
XXX
     48 	5.019874e+00
     49 	5.019850e+00
     50 	5.019666e+00
     51 	5.019626e+00
XXX
XXX
     52 	5.019556e+00
XXX
XXX
     53 	5.019529e+00
     54 	5.019505e+00
     55 	5.019502e+00
     56 	5.019501e+00
     57 	5.019500e+00
     58 	5.019381e+00
     59 	5.018639e+00
     60 	5.018618e+00
     61 	5.018404e+00
     62 	5.018242e+00
     63 	5.018046e+00
     64 	5.018046e+00
     65 	5.018046e+00
     66 	5.017979e+00
XXX
XXX
     67 	5.017645e+00
     68 	5.017217e+00
     69 	5.017132e+00
     70 	5.016711e+00
     71 	5.016438e+00
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     72 	5.016420e+00
     73 	5.016230e+00
     74 	5.015999e+00
     75 	5.015992e+00
     76 	5.015958e+00
     77 	5.015912e+00
XXX
XXX
XXX
XXX
XXX
     78 	5.015911e+00
XXX
XXX
XXX
XXX
XXX
XXX
     79 	5.015869e+00
     80 	5.015868e+00
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     81 	5.015836e+00
     82 	5.015804e+00
     83 	5.015741e+00
     84 	5.015712e+00
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     85 	5.015666e+00
     86 	5.015658e+00
     87 	5.015655e+00
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
XXX
     88 	5.015607e+00
XXX
XXX
XXX
XXX
     89 	5.015604e+00
XXX
XXX
XXX
XXX
XXX
     90 	5.015599e+00
     91 	5.015598e+00
     92 	5.015598e+00
     93 	5.015598e+00
     94 	5.015598e+00
     95 	5.015597e+00
     96 	5.015562e+00
XXX
XXX
     97 	5.015560e+00
XXX
     98 	5.015559e+00
     99 	5.015558e+00
    100 	5.015558e+00

Increasing Soft Maximum Generation Limit by 100 (MaxGenerations) to 200.
I'm doing this because the fitness is still improving.

    101 	5.015558e+00
    102 	5.015557e+00
    103 	5.015557e+00
    104 	5.015557e+00
    105 	5.015557e+00
    106 	5.015557e+00
    107 	5.015557e+00
    108 	5.015556e+00
    109 	5.015556e+00
    110 	5.015556e+00
    111 	5.015556e+00
    112 	5.015556e+00
    113 	5.015556e+00
    114 	5.015556e+00
    115 	5.015556e+00
    116 	5.015556e+00
    117 	5.015556e+00
    118 	5.015556e+00
    119 	5.015556e+00
    120 	5.015556e+00
    121 	5.015556e+00
    122 	5.015556e+00
    123 	5.015556e+00
    124 	5.015556e+00
    125 	5.015556e+00
    126 	5.015556e+00
    127 	5.015556e+00
    128 	5.015556e+00
    129 	5.015556e+00
    130 	5.015556e+00
    131 	5.015556e+00
    132 	5.015556e+00
    133 	5.015556e+00
    134 	5.015556e+00
    135 	5.015556e+00
    136 	5.015556e+00
    137 	5.015556e+00
    138 	5.015556e+00
    139 	5.015556e+00
    140 	5.015556e+00
    141 	5.015556e+00
    142 	5.015556e+00
    147 	5.015556e+00
    148 	5.015556e+00
    149 	5.015556e+00
    150 	5.015556e+00
    151 	5.015556e+00
    152 	5.015556e+00
    153 	5.015556e+00
    154 	5.015556e+00
    157 	5.015556e+00
    158 	5.015556e+00
    159 	5.015556e+00
XXX
XXX
    160 	5.015555e+00
    161 	5.015554e+00
XXX
    162 	5.014964e+00
    163 	5.014819e+00
XXX
XXX
    164 	5.014776e+00
    165 	5.014722e+00
    166 	5.014714e+00
    167 	5.014701e+00
    168 	5.014686e+00

WARNING: Turning Off MemoryMatrix because memory usage is too great.

    169 	5.014626e+00
XXX
XXX
XXX
XXX
    170 	5.014529e+00
    171 	5.014522e+00
    172 	5.014519e+00
XXX
XXX
XXX
    173 	5.014518e+00
    174 	5.014513e+00
    175 	5.014511e+00
    176 	5.014510e+00
    177 	5.014509e+00
    178 	5.014507e+00
    179 	5.014505e+00
    180 	5.014500e+00
    181 	5.014498e+00
    182 	5.014485e+00
    183 	5.014452e+00
    184 	5.014410e+00
    185 	5.014383e+00
    186 	5.014337e+00
    187 	5.014313e+00
    188 	5.014240e+00
    189 	5.014213e+00
    190 	5.014183e+00
    191 	5.014085e+00
    192 	5.014052e+00
    193 	5.013990e+00
    194 	5.013963e+00
    195 	5.013945e+00
    196 	5.013880e+00
    197 	5.013868e+00
    198 	5.013857e+00
    199 	5.013857e+00
    200 	5.013857e+00

Increasing Soft Maximum Generation Limit by 200 (MaxGenerations) to 400.
I'm doing this because the fitness is still improving.

    201 	5.013857e+00
    202 	5.013857e+00
    203 	5.013857e+00
    204 	5.013857e+00
    205 	5.013857e+00
    206 	5.013857e+00
    207 	5.013857e+00
    208 	5.013857e+00
    209 	5.013857e+00
    210 	5.013857e+00
    211 	5.013856e+00
    212 	5.013856e+00
    213 	5.013856e+00
    214 	5.013856e+00
    215 	5.013856e+00
    216 	5.013856e+00
    217 	5.013856e+00
    218 	5.013856e+00
    219 	5.013856e+00
    220 	5.013856e+00
    221 	5.013856e+00
    222 	5.013856e+00
    223 	5.013856e+00
    224 	5.013856e+00
    225 	5.013856e+00
    226 	5.013856e+00
    227 	5.013856e+00
    228 	5.013856e+00
    229 	5.013856e+00
    230 	5.013856e+00
    231 	5.013855e+00
    232 	5.013854e+00
    233 	5.013849e+00
    234 	5.013847e+00
    235 	5.013845e+00
    236 	5.013845e+00
    237 	5.013845e+00
    238 	5.013845e+00
    239 	5.013845e+00
    240 	5.013845e+00
    241 	5.013845e+00
    242 	5.013845e+00
    243 	5.013844e+00
    244 	5.013844e+00
    245 	5.013844e+00
    246 	5.013844e+00
    247 	5.013844e+00
    248 	5.013844e+00
    249 	5.013844e+00
    250 	5.013844e+00
    251 	5.013844e+00
    252 	5.013844e+00
    253 	5.013843e+00
    254 	5.013843e+00
    255 	5.013843e+00
    256 	5.013843e+00
    257 	5.013843e+00
    258 	5.013843e+00
    259 	5.013843e+00
    260 	5.013843e+00
    261 	5.013843e+00
    262 	5.013842e+00
    263 	5.013842e+00
    264 	5.013842e+00
    265 	5.013842e+00
    266 	5.013842e+00
    267 	5.013842e+00
    268 	5.013842e+00
    269 	5.013842e+00
    270 	5.013842e+00
    271 	5.013842e+00
    272 	5.013842e+00
    273 	5.013842e+00
    274 	5.013842e+00
    275 	5.013842e+00
    276 	5.013842e+00
    277 	5.013842e+00
    278 	5.013842e+00
    279 	5.013842e+00
    280 	5.013842e+00
    281 	5.013842e+00
    282 	5.013842e+00
    283 	5.013842e+00
    284 	5.013842e+00
    285 	5.013842e+00
    286 	5.013842e+00
    287 	5.013842e+00
    288 	5.013842e+00
    289 	5.013842e+00
    290 	5.013842e+00
    291 	5.013842e+00
    292 	5.013842e+00
    293 	5.013842e+00
    294 	5.013842e+00
    295 	5.013842e+00
    296 	5.013842e+00
    297 	5.013842e+00
    300 	5.013842e+00
    301 	5.013842e+00

Soft Generation Wait Limit Hit.
No Improvement in 30 Generations

Best Fit Found at Generation 301
Fit Value = 5.013842e+00


Parameters at the Solution (value, gradient):

 X[ 1] :	-5.963136e-01	G[ 1] :	-1.986686e+00
 X[ 2] :	4.627812e+00	G[ 2] :	6.480881e-01
 X[ 3] :	1.314116e+01	G[ 3] :	5.398500e-01
 X[ 4] :	-8.111620e-01	G[ 4] :	-1.015166e+00
 X[ 5] :	-3.941135e-02	G[ 5] :	-3.080531e+01
 X[ 6] :	7.225514e-01	G[ 6] :	-1.704932e+01
 X[ 7] :	1.281075e+00	G[ 7] :	-2.833139e+01
 X[ 8] :	6.387604e+02	G[ 8] :	-2.386708e-02
 X[ 9] :	-1.131333e+00	G[ 9] :	-2.909423e-01
 X[10] :	4.376772e-02	G[10] :	-1.237820e+02
 X[11] :	3.339252e+00	G[11] :	5.782831e+01
 X[12] :	3.057390e+00	G[12] :	1.994568e+01
 X[13] :	1.085533e+00	G[13] :	7.716009e+00
 X[14] :	1.457520e+00	G[14] :	-1.234986e+02
 X[15] :	2.474029e-02	G[15] :	1.025029e+02
 X[16] :	4.081959e+00	G[16] :	-4.501170e+01
 X[17] :	1.732660e+00	G[17] :	-2.927064e+01
 X[18] :	2.601479e+00	G[18] :	4.389799e+00
 X[19] :	2.762556e+00	G[19] :	9.191509e+01
 X[20] :	-7.821640e-03	G[20] :	1.619273e+03

Solution Found Generation 301
Number of Generations Run 301


Fri Mar  5 01:05:04 2004
Total run time : 0 hours 19 minutes and 17 seconds

(iter 2 ): LQD Results:
                                                            buchanan00      nader00     gore00      bush00 other00
(Intercept)/(Intercept)/(Intercept)/(Intercept)/NA         -0.59631359   0.72255141 3.33925188  4.08195934       0
pperot96/pdemnader96/pdem96/prep96/NA                       4.62781232   1.28107450 3.05739042  1.73265972       0
pChangeRepRV00/pGreenRV00/pChangeDemRV00/pChangeRepRV00/NA 13.14115625 638.76042224 1.08553317  2.60147905       0
pcuban00/pcuban00/pcuban00/pcuban00/NA                     -0.81116197  -1.13133337 1.45751985  2.76255633       0
pr1...1./pr1...2./pr1...3./pr1...4./NA                     -0.03941135   0.04376772 0.02474029 -0.00782164       0

(iter 2 ): LQD sigma: 5.013842 

WARNING: New fit is worse.  Stopping.
It may be a good idea to restart with a larger GENOUD population.
> 
> summary(muFLfull1, weights=TRUE)

Choice 1 : buchanan00 Estimates and SE:
                   Est SE.Sand SE.OPG SE.Hess t.val.Sand
(Intercept)    -0.6030   0.147 0.0581 0.03950      -4.11
pperot96        4.6700   1.310 0.4760 0.34100       3.55
pChangeRepRV00 13.2000   1.240 0.5720 0.36500      10.70
pcuban00       -1.0500   0.448 1.4000 0.20200      -2.35
pr1...1.       -0.0507   0.024 0.0140 0.00738      -2.11


Choice 2 : nader00 Estimates and SE:
                    Est SE.Sand  SE.OPG  SE.Hess t.val.Sand
(Intercept)      0.7250  0.1810 0.05770  0.03840       4.00
pdemnader96      1.2800  0.4500 0.13800  0.08680       2.86
pGreenRV00     638.0000 76.6000 0.03060 15.30000       8.33
pcuban00        -1.0300  0.2750 1.82000  0.13400      -3.73
pr1...2.         0.0359  0.0141 0.00524  0.00318       2.55


Choice 3 : gore00 Estimates and SE:
                 Est SE.Sand SE.OPG SE.Hess t.val.Sand
(Intercept)    3.350  0.1500 0.0294 0.02590      22.30
pdem96         3.050  0.3110 0.0684 0.05300       9.82
pChangeDemRV00 1.090  0.9140 0.1770 0.16000       1.19
pcuban00       1.470  0.2360 0.9080 0.11100       6.23
pr1...3.       0.022  0.0103 0.0026 0.00209       2.13


Choice 4 : bush00 Estimates and SE:
                    Est SE.Sand SE.OPG SE.Hess t.val.Sand
(Intercept)     4.08000 0.13800 0.0309 0.02340      29.50
prep96          1.73000 0.28600 0.0612 0.04970       6.05
pChangeRepRV00  2.61000 0.71000 0.1550 0.13100       3.68
pcuban00        2.76000 0.22100 0.9240 0.10900      12.50
pr1...4.       -0.00795 0.00743 0.0020 0.00163      -1.07


Choice 5 : other00 Estimates and SE:
               Est SE.Sand SE.OPG SE.Hess t.val.Sand
NA               0       0      0       0        NaN
NA               0       0      0       0        NaN
NA               0       0      0       0        NaN
NA               0       0      0       0        NaN
NA               0       0      0       0        NaN


LQD sigma: 4.966367 
TANH sigma: 4.318946 

Number of Observations: 67 
Number of observations with at least one zero weight: 8 
Number of zero weights: 11 
TANH: weights
   name weights:buchanan00 weights:nader00 weights:gore00 weights:bush00
1     1           1.000000        1.000000      1.0000000       0.000000
2     2           1.000000        1.000000      1.0000000       1.000000
3     3           1.000000        1.000000      1.0000000       1.000000
4     4           1.000000        1.000000      1.0000000       1.000000
5     5           1.000000        1.000000      0.8479780       1.000000
6     6           1.000000        0.000000      0.0000000       1.000000
7     7           0.812443        1.000000      1.0000000       1.000000
8     8           1.000000        1.000000      1.0000000       1.000000
9     9           1.000000        0.548565      1.0000000       1.000000
10   10           1.000000        1.000000      1.0000000       1.000000
11   11           1.000000        1.000000      1.0000000       1.000000
12   12           1.000000        1.000000      1.0000000       1.000000
13   13           1.000000        1.000000      1.0000000       1.000000
14   14           1.000000        1.000000      1.0000000       1.000000
15   15           1.000000        1.000000      1.0000000       1.000000
16   16           1.000000        0.000000      0.0000000       1.000000
17   17           1.000000        1.000000      0.2445280       1.000000
18   18           1.000000        1.000000      1.0000000       1.000000
19   19           1.000000        1.000000      1.0000000       1.000000
20   20           1.000000        1.000000      1.0000000       1.000000
21   21           1.000000        1.000000      1.0000000       1.000000
22   22           1.000000        1.000000      1.0000000       1.000000
23   23           1.000000        1.000000      1.0000000       1.000000
24   24           1.000000        1.000000      1.0000000       1.000000
25   25           1.000000        1.000000      1.0000000       1.000000
26   26           1.000000        1.000000      1.0000000       1.000000
27   27           1.000000        0.946213      1.0000000       1.000000
28   28           1.000000        1.000000      1.0000000       1.000000
29   29           1.000000        1.000000      0.9072720       1.000000
30   30           1.000000        1.000000      1.0000000       1.000000
31   31           1.000000        1.000000      1.0000000       1.000000
32   32           1.000000        1.000000      1.0000000       1.000000
33   33           1.000000        1.000000      1.0000000       1.000000
34   34           1.000000        1.000000      1.0000000       1.000000
35   35           1.000000        1.000000      1.0000000       1.000000
36   36           1.000000        1.000000      1.0000000       1.000000
37   37           1.000000        0.218009      0.0000000       1.000000
38   38           1.000000        1.000000      1.0000000       1.000000
39   39           1.000000        1.000000      1.0000000       1.000000
40   40           1.000000        1.000000      1.0000000       1.000000
41   41           1.000000        1.000000      1.0000000       1.000000
42   42           1.000000        1.000000      0.1896080       0.262881
43   43           1.000000        1.000000      0.0403284       1.000000
44   44           1.000000        0.213169      1.0000000       1.000000
45   45           1.000000        1.000000      1.0000000       1.000000
46   46           1.000000        1.000000      1.0000000       1.000000
47   47           1.000000        1.000000      1.0000000       1.000000
48   48           1.000000        0.000000      1.0000000       1.000000
49   49           1.000000        0.926798      1.0000000       0.922739
50   50           0.000000        0.000000      1.0000000       1.000000
51   51           1.000000        0.445467      0.0000000       1.000000
52   52           1.000000        1.000000      0.0000000       1.000000
53   53           1.000000        0.411884      1.0000000       1.000000
54   54           1.000000        1.000000      1.0000000       1.000000
55   55           1.000000        1.000000      1.0000000       1.000000
56   56           1.000000        1.000000      0.7648670       1.000000
57   57           1.000000        1.000000      0.4807590       1.000000
58   58           1.000000        0.958557      1.0000000       1.000000
59   59           1.000000        0.495407      1.0000000       1.000000
60   60           1.000000        1.000000      1.0000000       1.000000
61   61           1.000000        1.000000      1.0000000       1.000000
62   62           1.000000        1.000000      1.0000000       1.000000
63   63           1.000000        1.000000      1.0000000       1.000000
64   64           1.000000        1.000000      0.6927510       1.000000
65   65           1.000000        1.000000      1.0000000       1.000000
66   66           1.000000        1.000000      1.0000000       1.000000
67   67           1.000000        1.000000      1.0000000       1.000000
> 
> save(muFLfull1, file="FLfull1.a.RData");
> proc.time()
[1] 2032.28  104.59 2148.20    0.00    0.01
> 
